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This paper presents a method to calculate the dynamic coefficients of seals with 
arbitrary geometry. To describe the turbulent flow the Navier-Stokes equations are 
used in conjunction with the k-e turbulence model. These equations are solved by a 
full 3-dimensional finite-difference procedure instead of the normally used 
perturbation analysis. The time dependence of the equations is introduced by 
working with a coordinate system rotating with the precession frequency of the 
shaft. The results of this theory are compared with coefficients calculated by a 
perturbation analysis and with experimental results. 


INTRODUCTION 

During the last years it has become evident that it is important to include the 
fluid forces caused by seals when predicting the dynamic behavior of turbopumps. To 
calculate these forces and the dynamic coefficients which are normally used to 
describe them (eq. 1) 
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several methods which are based either on the so-called "bulk-flow” theories /1,2/ 
or directly on the Navier-Stokes equations /3.4/ have been published. A common 
feature of these methods, developed for straight pump seals /l/ straight gas seals 
/2,4/ or grooved seals /3/, is that they are all based on a perturbation analysis 
to determine the dynamic coefficients. But the perturbation analysis requires many 
assumptions, e.g. 

1. It is assumed that the shaft moves on small orbits around the centric position, 
so that a perturbation analysis can be used for all flow variables 

2 

$ = ♦ o + e ^i + e ^2 

and all terms with power of e greater than one can be neglected in the 
equations without loss of accuracy (e = perturbation parameter). 


211 



2. The change of the perturbation flow variables in the circumferential direction 
can be described by sine and cosine functions. 

*1 = * lc cos0 + * ls sin0 

3. The change in time can be described by 


♦l = *1 6 


iQt 


because the shaft moves on a circular orbit. 

To check how these assumptions effect the results we have developed a 3-dimensional 
finite-difference procedure to calculate the dynamic coefficients. The only 
assumptions in this theory are that the turbulence can be described by the k-e 
turbulence model and that the shaft moves on circular orbits around the seal 
center . 


GOVERNING EQUATIONS 

To describe the turbulent flow in a seal we have the Navier-Stokes equations and 
the continuity equation. The turbulent stresses occurring in the fluid can be 
handled like laminar stresses by introducing a turbulent viscosity. The turbulent 
and the laminar viscosity are then summed up to an effective viscosity p & 


^e = **1 + ^t 


To describe p^ the k-e turbulence model /5,6/ is used because it is simple and 
often used to calculate the turbulent flow in seals /7.8.9.10/. 


„ k 

P .SB C P 

*t p e 


All these equations can be represented in the following form 
3* d , ^ d fr 3*> . 13, .v 13, r 3*. A 13, 1 3, r ia *\ o 

PaF + d^ pu<p) ■ d^ r * dJ + 73F (rpv * ) ■ 73F< rr * aF* + Fae^*" F30< r * Fa0 } = s * 
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Table 1: Source terms of the Navier-Stokes equations, the 
continuity equation and the equations of the k-e 
model . 

(constants of k-e model are given in Appendix A) 


To determine the dynamic coefficients we assume that the shaft moves on a circular 
orbit with precession frequency 0 around the seal center. Since this would normally 
result in a time dependent problem we introduce a rotating coordinate system which 
is fixed at the shaft center (Fig. 1). In this system the flow is stationary. 

Due to the rotating coordinate system centrifugal- and coriolis-forces occur in the 
equations for the radial and circumferential momentum /ll/, which are taken into 
consideration by a modification of the source terms. 

S 1 = S + 0^ + 20w 

V V 

s' = S -20 v 
w w 

So the final form of the equations is given if S^ and S^ in Table 1 are replaced by 

S' and S 1 . 
v w 
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Fig. 1: Rotating coordinate system 


BOUNDARY CONDITIONS 

The above given equations are solved in conjunction with the following boundary 
conditions (Fig. 2) 



Fig. 2 : Locations where the boundary-conditions must be 
specified 
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Entrance: p B(0) + | pu^ (0) (l+f) = p A 

^ : P C (0) = 0 

Stator : u=0 v = 0 w = -fir 

Rotor : u = 0 v = 0 w =(u-fi)r 

u D is the average axial entrance velocity specified for every grid-plane. 

15 

For k and e the standard conditions of the k-e model /5,6/ are used at the walls. 


FINITE-DIFFERENCE-PROCEDURE 

This system of equations with the corresponding boundary conditions is solved by a 
3-dimensional finite-difference procedure, based on the method published by Gosman 
and Pun /12/. The seal is discretized by a grid (Fig. 3) and the variables are 
calculated at the nodes. To determine the pressure we use the PISO /13/ algorithm 
instead of SIMPLE /14/. 



Fig. 3: Mesh arrangement in the seal 


DYNAMIC COEFFICIENTS 

As result of our solution procedure we get the pressure distribution in the seal 
and by a pressure integration in axial and circumferential direction, the forces. To 
simplify the integration we consider the case when the rotating coordinate system 

y, z coincides with the stationary system y, z. 
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When calculating the forces for 3 precession frequencies Q = Ou. Q = lu and fi = 2io 
we can determine the dynamic coefficients. To save computation time we calculate k, 
e and only for fi = 0 and keep it constant then for 0 = lw and 0 = 2 o>. 


RESULTS 


1. Example 


First we compare the results of the 3-dimensional theory with those of a method, 
based on a perturbation analysis in the Navier-Stokes equations and in addition 
with experimental values. 

For a straight seal with the following data the results are shown in Fig. 5. 


L = 23.5 ran 

Tj — 23.5 ran 

C = 0.2 ran 

o 

w(0, 6,2000 RPM)/r i <j 
w(0, 6.4000 RPMJ/r^ 
w(0, 6,6000 RPM)/r i <o 


0.13 

0.17 

0.19 


Hj = 0.7*10 2 3 Ns/m 3 
P = 996.0 kg/m 3 

u = 16.46 m/s 

f (2000 RPM) = 0.35 
f (4000 RPM) = 0.37 
f (6000 RPM) =0.38 


2. Example 


For a grooved seal (Fig. 4) 
groove depth 1^ (Fig. 6.7). 


the dynamic coefficients are shown as a function of the 



Fig. 4: Geometry of grooved seals 
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The seal data are 


L = 23.5 

mm 

= 0.7*10 3 Ns/m 3 

r. = 23.5 

mm 

p = 996.0 kg/m 3 

C = 0.2 
o 

mm 


w(0,9.u = 

14.11 m/s )/t ± u = 0.20 

f (14.11 m/s) = 0.38 

w(O,0,u = 

11.76 m/sJ/r.G) = 0.22 

f (11.76 m/s) = 0.39 


In Fig. (8,9) the grid, the axial and radial velocity, the circumferential velocity 
and the pressure are shown for a groove depth of 0.5 mm and an average axial 

velocity of u = 14.11 m/s. 


3. Example 

In some further calculations we have investigated the influence of grooves on the 
rotor and stator for the seal shown in Fig. 10. 

The seal data are 


L = 35.0 

mm 

**1 

= 0.7*10 3 Ns/m 3 

r j = 23.7 

mm 

p 

= 996 . 0 kg/m 3 

C = 0.2 
o 

mm 

n 

= 4000 RPM 

£ = 0.5 


w(O,0) 

= 0. 5*r j<d 


groove depth on rotor and stator : 0.4 mm 

total pressure loss : 0. 8 Npa 

radius of shaft orbit : r = C /40 

o o 

In Fig. 10 the total stiffness coefficients of the seal, and the portions developed 
in each part of the seal are shown. 
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Fig. 5: Comparison of the dynamic coefficients for a 
straight seal. 
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Fig. 6: Comparison of the dynamic coefficients for grooved 
seals. 
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Fig. 7: Comparison of the dynamic coefficients for grooved 






Fig. 8: Grid and axial and radial velocity for a grooved 
seal 
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Fig. 9= Pressure distribution and circumferential velocity 
for a grooved seal. 
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100000 




Fig. 10: Stiffness coefficients for a seal with grooves on 
rotor and stator 

From that diagram we can draw the following interesting conclusions: 

1. Although the pressure loss for a centric shaft position has the same magnitude 
for every land part, these lands develop different contributions to the total 
coefficients. 

2. Although the clearance in the chambers is 5-times greater than in the land 
parts, the forces in the chambers can’t be neglected. 

3. The chambers have a strong destabilizing effect, because they cause positive 
radial forces and big positive tangential forces, (positive forces have the 
direction of the z-y axes in Fig. 1) 


4. Example 


We made further test calculations for the seal arrangement shown in Fig. 11. 
The seal data for this example are 


Pj = 0.7e 3 Ns/m^ 
p = 1000 kg/m 3 

total pressure loss : 0.347 Mpa 

Radius of shaft orbit : r = C /40 

o o 


f = 0.5 

w(0,e)/r i( j = 0.5 
n = 2000 RPM 
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The result is surprising, because seal 1 yields a negative direct stiffness instead 
of the expected positive. 

This can be explained if one looks at the pressure loss in the seal arrangement. 




D 



Fig. 13: Pressure loss for 

centric shaft 
position 


Pressure loss in the 
plane with nearest 
gap in seal 1 


Pressure loss in the 
plane with widest 
gap in seal 1 


The main pressure loss is caused by the entrance loss of the flow from the chamber 
into seal 2 (Fig. 13). If now, the rotating part moves in the redirection the increase 
of the clearance of seal 2 (in the plane considered) results in a sharp drop in the 
pressure loss B-C. And this drop in the pressure loss and the rise in the opposite 
plane is responsible for the negative value of k in seal 1. 


CONCLUSIONS 

The first two examples show that the perturbation analysis yields good results in 
comparison with the 3-dimensional theory although it requires only a fraction of 
the calculation time and storage needed for that method. On the other hand it is 
only possible to determine the coefficients of look- through seals with the 
perturbation analysis, while there are no restrictions concerning the geometry for 
the 3-dimensional procedure. 

Example 3 and 4 clearly demonstrate that in an arrangement of several seals, the 
separated consideration of the single seals, even with big chambers between them, 
may lead to totally wrong results. 
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NOMENCLATURE 


F , F 
z y 

K. k 

D, d 

M» m 

u,v,w 

P 

k 

e 

■V^t 

p 

t 

x.r.0 

G 

a k ,a e 

V c r c 2 


e = r /C 
o o 

w 

0 

S 

L 

r. 


Forces on the shaft in z and y direction 

direct and cross-coupling stiffness 
direct and cross-coupled damping 
direct and cross-coupling inertia 
axial, radial and circumferential velocity 
pressure 

turbulence energy 
energy dissipation 

effective, laminar and turbulent viscosity 

density 

time 

axial, radial and circumferential coordinate 
production term in k-e-model 
Constants of the k-e-model 

Constants of the k-e-model 

general variable standing for u.v.w.p.k e 
general source term 

seal clearance 

radius of the precession motion of the shaft 

perturbation parameter 

rotational frequency of the shaft 
precession frequency of the shaft 
entrance lost-coefficient 
Length of the seal 
radius of the rotor (shaft) 

radius of the stator 


APPENDIX A: Constants and production term of the k-e model 


C^ = 0,09 Cj = 1,44 C 2 = 1,92 = 1,0 o fe = 1,3 


G = u 


o ,, 9 ** 2 ,du 2 I aw A v 2 . fdv du 2 . A dv . dw _ w 2 . ,dw . 1 du. J 

^dr^ ^dx^ r d 0 r^ ^dx dr^ ^r d 0 dr r^ ^dx r § 0 ^ 
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